Computational Power and Correlation in Quantum Computational Tensor Network 
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We investigate relationships between computational power and correlation in resource states for quantum 
computational tensor network, which is a general framework for measurement-based quantum computation. We 
find that if the size of resource states is finite, not all resource states allow correct projective measurements in 
the correlation space, which is related to non-vanishing two-point correlations in the resource states. On the 
other hand, for infinite-size resource states, we can always implement correct projective measurements if the 
resource state can simulate arbitrary single-qubit rotations, since such a resource state exhibits exponentially 
decaying two-point correlations. This implies that a many-body state whose two-point correlation cannot be 
upperbounded by an exponentially decaying function cannot simulate arbitrary single-qubit rotations. 

PACS numbers: 03.67 .Lx 



I. INTRODUCTION 

Quantum many-body states, which have long been the cen- 
tral research objects in condensed matter physics and sta- 
tistical physics, are now attracting the renewed interest in 
quantum information science. In particular, in measurement- 
based quantum computation (MBQC) Illy], the role of quan- 
tum many-body states played in quantum computation is very 
clear: once the special resource state, which is called the clus- 
ter state, is prepared, universal quantum computation is pos- 
sible with adaptive local measurements on each qubit. This 
clear separation between the resource preparation and the exe- 
cution of computation itself has also relaxed the requirements 
for experimental realization of quantum computation fl-Qj]]. 
Recently, the concept of quantum computational tensor net- 
work (QCTN), which is a novel framework of MBQC on 
general many-body states, has been proposed lfl2l - [l4ll . The 
most innovative feature of this framework is that the resource 
state is represented in the matrix-product (or tensor-product) 
form il5l4l7ll and universal quantum computation is simulated 
in the virtual linear space, which is called the correlation space 
(CS) cull, where matrices (or tensors) live. 

How does the behavior of a physical quantity in a many- 
body resource state affect the computational power of a quan- 
tum computer? If one wants to construct a bridge between 
quantum many-body physics and quantum information sci- 
ence, one cannot avoid tackling such a question. For exam- 
ple, in Refs. (H H, the beautiful result was obtained that 
if a many-body resource state has too much entanglement, it 
is useless for MBQC. In Ref. [20], a criterion of the amount 
of entanglement for universal preparators in graph states was 
obtained. In Ref. l2lll . a relationship between the amount of 
entanglement and the gate fidelity in the one-way model was 
derived. These results have clarified important relationships 
between entanglement in many-body states and the computa- 
tional power of MBQC. 

Surprisingly, less attention has been paid to relationships 
between more traditional physical quantities, namely, a two- 



point correlation and the computational power of a quan- 
tum computer. Obviously, a two-point correlation is one of 
the most important research subjects in traditional quantum 
many-body physics, such as condensed matter physics and 
statistical physics. If we consider applications of quantum 
many-body states in condensed matter physics to quantum 
computation B22I1 . it is natural to ask how the behavior of 
a two-point correlation affects the computational power of a 
quantum computer. It has long been known in quantum in- 
formation science that two-point correlations exactly vanish 
within a finite distance in the cluster state. Such a correla- 
tion property prohibits the cluster state from being the exact 
ground state of any physical Hamiltonian Il23l425il . Recently, 
it was shown that the ground state of the Affleck, Kennedy, 
Lieb, and Tasaki (AKLT) model M26I1 . which exhibits the ex- 
ponentially decaying two-point correlation, enables univer- 
sal MBQC llT2l - [l4ll27ll . This result suggests that we might 
be able to further explore other resource states in the pool 
of quantum many-body states including ones which exhibit 
polynomially decaying or non decaying two-point correla- 
tions. Considering the fact that there are plenty of many-body 
states which exhibit such a long-range two-point correlation 
in condensed matter physics, it is very interesting to study 
whether we can use long-range-correlated many-body states 
as resources of MBQC. 

In this paper, we investigate relationships between compu- 
tational power and two-point correlations in QCTN consider- 
ing both finite and infinite-size resource states. Specifically, 
we here address under what conditions we can simulate pro- 
jective measurements and arbitrary single-qubit rotations in 
QCTN, both of which are important properties of a universal 
resource for MBQC j2§l . 

We find that projective measurements are not simulated cor- 
rectly in QCTN on certain resource states of finite size, such 
as the AKLT states with a specific boundary condition. We 
therefore introduce a class of finite-size matrix-product states 
(MPSs), namely, normalizable resources, which allow correct 
projective measurements in the CSs. Since two-point corre- 
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lations exactly vanish within a finite distance in normalizable 
resources, they are considered as cluster like resources. In 
order to perform correct projective measurements on finite re- 
source states with non vanishing two-point correlations, we 
introduce another class of MPSs. We also mention an alterna- 
tive way of simulating projective measurements indirectly by 
using the downloading method j29ll . 

We further show that if a resource state can simulate arbi- 
trary single-qubit rotations, it can also perform correct projec- 
tive measurements in the infinite-size limit, that is, it is asymp- 
totically normalizable. Accordingly, normalizability (at least 
asymptotic sense) is a necessary and sufficient condition for 
a universal quantum computational wire, where both correct 
projective measurements and arbitrary single-qubit rotations 
can be simulated. This result is obtained by showing a the- 
orem about convergence of a stochastic unitary map. Fur- 
thermore, this theorem also tells us that a many-body state 
whose two-point correlation cannot be upperbounded by an 
exponentially decaying function, such as the Greenberger- 
Horne-Zeilinger (GHZ), the Ising antiferromagnetic, and the 
W states, cannot be used to perform an arbitrary single-qubit 
rotation. These findings highlight interesting relationships be- 
tween computational power and correlation in many-body en- 
tangled states. 

This paper is organized as follows. In Sec. I, we briefly 
review QCTN IIL?l - [l4ll . which is the most general framework 
for MBQC. In Sec. II, we investigate simulations of projec- 
tive measurements in the QCTN framework and specify un- 
der what situation one can accurately simulate projective mea- 
surements in the resource states of finite size. In Sec. Ill, we 
extend the notion of normalizability to the resource states of 
infinite size, namely asymptotic normalizability, and address 
a fundamental relationship between asymptotic normalizabil- 
ity and computational power by providing a theorem which 
makes a bridge between the ability to perform an arbitrary 
single-qubit rotation and behavior of the two-point correla- 
tions. Section IV is devoted to the conclusion. 



II. QUANTUM COMPUTATIONAL TENSOR NETWORK 

Let us consider an MPS SCI]: 

\9(R,L))? = (R\A[i n ]---A[i 1 }\L)\i n ---i 1 ), 
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FIG. 1: (a) Matrix product state, (b) By measuring from the 1st to 
(k — l)th qudits, a single qubit rotation is simulated in the correla- 
tion space, \L) —¥ \ip). In the normalizable resource with a finite 
I, projective measurements are simulated correctly at every position 
which is more than I qudits away from the right boundary. 



by performing a projective measurement in the basis {\otj) = 



EfLo a J»l*)} = °>' 



1 and (a.i\a.j) — r5y) on a 



physical qudit. The postmeasurement state with the measure- 
ment outcome a, reads 



(R\A[i n ]---A[i 2 ]A[a j ]\L)\ 
= |*(i2,Z/))a®>/5Jl«j>. 



In ■■ •*2> <8> \OLj 



(1) 



where A[a s ] = £to a^Aft and \V) = A[atj]/ y/aj\L) 
with a,jl = A*[aj]A[aj]. The postmeasurement state can be 
understood that the boundary vector | L) is transformed to \L') 
by the single-qubit operation A[ctj}/ y/aj. So far, the normal- 
ization of an MPS has not been considered properly. In order 
to address this point, we define the normalization factor as 

f n (\R),\L))^(R\(A° n \L)(L\)\R), 

where Ap = YliZn A[i]pA^[i\. We denote the normalized 
MPS as 

\<z(R,m = MR,m/\/u\R),\L)). 



where A[i k ] and \R), \L) ((R\R)=(L\L) = 1) are a D x D 
matrix and /^-dimensional boundary vectors IIL2T - [l4ll . respec- 
tively, which live in the virtual Hilbert space, the so-called 
correlation space (CS) IHG1 [see Fig. Q](a)]. The \i k ) 
(ik = 0, • • • d — 1 and = 5ij) denotes the fcth physi- 
cal qudit. Since we are interested in the simulation of unitary 
operations in the CS, we assume that the MPS matrices A[i] 
of a resource state can be described as = cj with 

positive real numbers Cj by choosing the physical basis prop- 

erly, where Eto c * = 1 ( Le - £to A^[i]A[i] = I). We 
here consider the case of D = 2, that is, the CS simulates 
a single qubit. We can simulate a quantum gate in the CS 



HI. SIMULATIONS OF PROJECTIVE MEASUREMENTS 

A. Condition for correct projective measurements 

Let us consider a situation where from the first to (k — l)th 
physical qudits are measured in order to implement a single- 
qubit rotation in the CS, and the boundary vector \L) is trans- 
formed to [see Fig. |T](b)]: 



(R\A[i n ] • • • A[i k ]\il>)\i n • • • i k ) 

y/f n -k+l(\R),m 
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Then, we want to perform a projective measurement on the 
state \tp) in the CS in a basis {|0), |1)}. To this end, without 
loss of generality, the MPS matrix A[i] can be written as 

A[0] = v / W2(l<M<o| + |<M(i|) 

A[l] = yfej2(\fo)(0\ - \<t>i)(l\) 

by choosing a proper physical basis, where |</>o) and \<f>\) 
are orthonormal states and c m is a real number. If we do 
the measurement M = {|mo) = (|0) + |1))/V2, |mi) = 
(|0) - \1))/V2, |2), • • • , \d - 1)} on the fcth physical qudit 
and obtain the outcome rrij (j = 0, 1), the probability pj of 
obtaining rrij (j = 0, 1) is calculated to be 

_ c m J n _ k {\R),\4>m3\^)\ 2 

Pi f n . k+1 (\R),m ■ 

If the above operation properly simulates the projective mea- 
surement in the CS, Po/pi = |(0|V>)| 2 /I(1| 1 /')| 2 must be satis- 
fied. This leads to 

f n -k(\R),\<h))=f n -k(\R),\4>i)). (2) 

If a resource state does not satisfy condition @, correct pro- 
jective measurements are not simulated in the CS. 

B. Normalizable resources 

In order to perform correct QCTN, we have to find MPSs 
which satisfy condition (0. Let us first consider a sufficient 
condition of Eq. (0. We call a resource state normalizable, if 
there exists a natural number I (I < n) such that 

/ n _ fc (|-R),|V)) = i/2, 

for all \ip), \R) and all fc < n — I. Obviously normalizable 
resources satisfy condition (O as long as n — k > I, which 
means that we can perform correct projective measurements 
in the CS at every position which is more than I qudits away 
from the left boundary. Later we will see that normalizability 
in the limit of I — > oo is also a necessary condition for an 
arbitrary single-qubit rotation. 

A resource state is normalizable, if and only if the map A° l 
(I < n — k) is a completely depolarizing channel: 

3 

£p = l/i^Uo-tVpV^a.^, 

i=0 

where at are the Pauli matrices (<tq = I), and U, V are ar- 
bitrary unitary operators (see Appendix A). If a finite-size re- 
source state (n < oo) is normalizable, I is also finite. Let 
us first consider the case of I = 1. In this case, A = £, 
which means that the dimension of a physical system is at 
least d = 4, and in such a case the MPS matrices are given 
by A[i] = UaiV/2. It is obvious that such a resource state 
with properly chosen U and V allows an arbitrary single- 
qubit rotation. Next we consider the case of I = 2. The 
minimum dimension of a qudit which satisfies A 02 = £ is 



d = 2, where the MPS matrices are given, for example, by 
v4[i]j4[j] = U(Ji + 2jV = 0, 1). Specifically, if we choose 
U and V in such a way that UaiW = Ug%U' = —i<J2, 
Ua-zU^ = cr 3 and V = U\ the MPS matrices become 
A[0] = H/V2 and = HZ/V2, which is the matrices 
for the cluster state U 12, [T^-H • Since the above two- and 
four-dimensional resource states are normalizable with a finite 
I, we can simulate correct projective measurements in the CSs 
even if their size is finite. 

The above definition of normalizability also has an impor- 
tant physical meaning: if a resource state is normalizable 
with a finite I (i.e., A° l = £), any two-point correlation 
(°fe°fc+r)* " (0£)*(Otk->* vanishes for all r > Z (see 
Appendix B), where (9^' b is an arbitrary operator on the fcth 
qudit, and (O)* = (*(L, R)\0\V{L, R))^. Thus the pos- 
sibility of projective measurements in QCTN with finite-size 
resource states is related to their vanishing two-point correla- 
tions. In the sense that two-point correlations exactly vanish 
with a finite distance, the normalizable resources with a finite 
I can be viewed as cluster like resource states 01 Q]. Then, 
it is natural to ask, can we simulate correct projective mea- 
surements on the resource states with non vanishing two-point 
correlations? 



C. A class for finite QCTN 

In order to answer that question, we here find a class of 
finite-size MPSs, which is not necessarily normalizable with 
a finite I, but can satisfy condition (0 with a finite / by prop- 
erly choosing the boundary condition \ R). This class of MPSs 
includes a wide variety of resource states, such as the AKLT 
states with specific boundary conditions |[L?l - [T4ll27l[3lll . As a 
canonical form of MPSs in the class, we take the matrices as 

A[0] - ^j2e l0 « Y ,A[l} = ^j2e^ Y Z, 

for j = 2, • • • , d - 1, where Bj e {I, X, Y, Z} and c m + 

Sj=2 c j = 1' i- e -' X)j=o ^H^H = !■ We can easily show 
that this class of MPSs with properly chosen 8i is capable of 
performing an arbitrary single-qubit rotation. Since A[i] con- 
sists of the Pauli matrices and rotations about the Y axis, the 
map A transforms the Pauli matrices as follows 

AY = J22Y, AX = j n X + ji 3 Z, AZ = l31 X + l33 Z, 

where 7^ are some real numbers. Using this fact and 
\4>j) (<f>3 1 = I / 2 +(- 1 Y ( cos 26> Z-sin 28 X) /2, we can con- 
clude that the normalization factor takes the following form: 

f n - k QR),\(f> j )) = (R\(I/2 + glZ + g(X)\R), 

with #3 1 being some real numbers, which depend on n and 
fc in general. Then, if we choose \R)(R\ = (/ ± Y)/2, 
fn-k(\R),\4>o)) - f n - k (\R),\<h.)) = 1/2 for all fc < n. 
That is, by choosing the boundary vector as \R) = (I±Y)/2, 
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we can perform correct projective measurements at every po- 
sition [32]. Instead of Y and Z, we can also use other com- 
binations of two Pauli operators in the canonical form (|3j. 
Specifically, if we set d = 3, O = -tt/2, 6 2 = 0, B 2 = Z, 
and c m /2 = c 2 — 1/3, the matrices become A[0] = XZ/y/3, 
Ml] = X/y/3, and A[3] = Z/y/3, which is the AKLT state 
rf27ll (we can also show that the AKLT-like state lfT2T - [l4Tl is an 
MPS of this class). 



D. Download as an indirect way to perform projective 
measurements 



IV. ASYMPTOTIC NORMALIZABILITY AND 
COMPUTATIONAL POWER 



Next we investigate a relationship between normalizability 
and possibility of an arbitrary single-qubit rotation by extend- 
ing the notion of normalizability to the infinite-size resource 
states. We call a resource state asymptotically normalizable 
if it is normalizable only in the limit I — > oo. For example, 
in the case of the AKLT state [27], the matrices are given by 
A[0] = X/\/3, A[l] = XZ/V3, and A[2] = Z/\/3, and the 
normalization factor can be calculated to be 

3 

f n - k (\R), |V» = 1/2 + (-l/3)"- fc J2 ei<i2k,| J2)/2, 



We here explore another way to simulate projective mea- 
surements indirectly. In Ref. [29], a method has been pro- 
posed to download a quantum state from the CS to the phys- 
ical space. In their method, we first project the fcth qudit of 
\^(R, to the subspace spanned by {|mo), (if we 

fail the projection, we repeat it again until it succeeds). Next 
we rotate the state \fa) to \rrij) (j = 0, 1) via MBQC. Then 
we now have 



^ \9(R, mj ))Z® 

3=0,1 



' c u c m f n -k>+i(\R), \rrij)) 
f n - k+1 {\R)M)) 



(j'IV')Itoj) 



where k' > k + 1, and the factor c u comes from the mea- 



surements for the rotation | faj 



in the CS. Finall y w e 



perform the projective measurement M. on the fc'th qudit 113011 . 
The postmeasurement state corresponding to the measurement 
outcome m; (I = 0, 1) is given by 

£ |*(i^j)>& + i®N) 
i=o,i 



,rryl CuC 2 m } n -k'{\R), \4>l)) 

V 2/„_fc + i(i?), V}) 



The coefficient of does not depend on j, and hence 

we can download the state from the CS to the physical space, 
even if the resource state is of finite size. This means that 
MBQC on finite-size resource states is universal state prepara- 
tor 112011 . Surprisingly, the above fact means that for finite-size 
resource states universal state preparator is more ubiquitous 
property than the simulator of a classical output. Of course we 
can perform projective measurements indirectly on a quantum 
state in the CS by downloading it to the physical space and 
measuring the downloaded qudit. However, as discussed in 
Ref. 12911 . a faithful downloading process takes a time step 
in proportion to the correlation length of the resource state, 
which results in additional operations. Since such additional 
operations can cause additional decoherence effects, it would 
be better to avoid doing download if we need only classical 
outputs. 



where \?p)(ip\ = 1/2 + Si=i e i <T i/2- I n me limit I — > oo, 
f n -k(\R), \ — > 1/2, and hence the AKLT state is asymp- 
totically normalizable. 

To make a bridge between asymptotic normalizability and 
possibility of an arbitrary single-qubit rotation, in the follow- 
ing, we show a theorem about convergence of completely- 
positive-trace-preserving (CPTP) maps, which allows a clas- 
sification of the MPSs of D = 2 with respect to their behavior 
of two-point correlations. 

Theorem. Suppose A is a stochastic unitary (hence CPTP) 
map. If A does not converge to £ in the limit of / — > oo 
(i.e., lirm^oo A ^ £), the Kraus operators of A can al- 
ways be described by choosing a proper basis as A[i]/ \fci = 
Xi 1 -^)/^Z{4> i ), where u 3 = ±1 and Z(<&) = e -^ z / 2 
with a certain angle fa . 

Proof. A stochastic unitary CPTP map A can be written in 
terms of the {^cri}^ =Q basis as a 4 x 4 matrix, 



A = 



/ 1 \ 


A 



where A is a 3 x 3 real matrix with 



< 1, where 



Oil no = max 



\Oi 



is the oo norm (see Appendix C for the matrix representation 
of the map ^4). 

If | \A\ |oo < 1> A ~ > £ m me limit of I — > oo, since 

Halloo <\\a\\L^o- 

Hence, in this case, the corresponding MPS state is at least 
asymptotically normalizable. 

On the other hand if ||^4.||oo = 1> there exist two vectors v 
and v' such that Av = v' and ||i>|| 2 = ||t»'|| 2 = v. We can 
always choose the vector v and v' as real vectors, because if 
v is a complex vector, the complex conjugate v* also subjects 
to Av* = 



Thus we can redefine a real vector v by v 



v*, which also satisfy Av 
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Let us consider two Hermitian operators 

1 3 

V = -=Y, V i a i> 

V = 

By their definitions, Tr[V] = Tr[V] = and Tr[V 2 ] = 
Tr[V' 2 ] = 2, and therefore both p„ = (7 + V)/2 and 
/V = (7 + V')/2 are pure states. Since = t/ (i.e., 
.4 V = V), Ap v = p v i ; that is, a pure state p v is mapped to 
another pure state p v > without any leakage (trace-preserving). 
On the other hand, we have assumed that the Kraus operator 
A[i] of the map A is proportional to a unitary operator. Thus 
the Kraus operator (MPS matrix) has to be written as 

-^=A[%\ = \i/)(v\+e iPi \v' ± )(v ± \ 



Ci 

with a certain angle /3j for all i, where \v)(v\ = p v and 

= p V ', and Iv^) and \v' ) are their orthogonals re- 
spectively. 

By redefining Pauli operators such that \v) and \v ± ) are the 
eigenstates of 03, the MPS matrix can be reformulated as 

4=A[i] = UZ(J3i), 

where Z[fii) = e -^ z / 2 and U = \v')(v\ + ^'^(v^l 
Then, the map A can be decomposed into U o Z, where 
Z = E i ^(A)p^ t (ft) c » and = Upf/t. The 3 x 3 ma- 
trices Z and W are also defined for Z and respectively. 
Specifically, the 3 x 3 matrix Z is written in the present basis 
as 



Z = 



where p = ^2 i cos(f3i)c i and q = Ej sm (ft) c i- By redefin- 
ing the (72 and (73 basis, the above matrix can be rewritten as 

/ Ai 
Z = A 2 
\ 1 

Note that 

l-p 2 -q 2 = ^c l (cos 2 ft + sin 2 A)-^c 4 (p 2 + (? 2 ) 

i i 

= ^c,[(cos/3 4 -p) 2 + (sin^ - q) 2 ] 

i 

> 0, 

where we have used cos 2 + sin 2 = 1, J2i c i = 1> and 
7^ j8j . This means that 




|Ai| = |A 2 | = V / F+9 5 <1- 



(The reason why we have assumed ^ /3j is that if j3\ = /3 2 , 
for example, then by changing the physical basis as 



|1> 

|2> 



V c l + c 2 
^2" 



|1> 



1) 



/ci 



\/cl + C 2 V C 1 + C 2 



|2>, 
|2>, 



we can make A[2] =0 in the MPS. ) 

The 3x3 matrix U, on the other hand, is an 0(3) rotation, 
since U is a unitary operation (i.e. it does not change ||v||). 
Let us consider a Hermitian matrix 

(ZUZ'^ZUZ' 
= Z'U T Z' 2 UZ> 


AZ r + Z'U T [ _ 
1-A 



uz' 



^.2 



U\ U 2 U\ U3U1 



XZ' + (1 - A)Z' u 2 ux u\ u 3 u 2 Z' 



U1U3 U 3 U 2 M3 



where m = U 3l , A = |Ai, 2 | 2 and Z' = diag(V% Vx, 1). The 
eigenvalues of the above matrix are given by A 2 and 



1 

2 

Since 



|(1 - A) 2 + 2A ± (1 - A) yj(l- X) 2 uj + AXu 



1 = U 2 + 1*2 + u\ > It 2 , 



r±. 



\r±\<~ 2 



(1- A) 2 + 2A + 



If |r+| = 1 or |r_ 
This implies that 



2A+(1-A)^(1- A) 2 +4A 
1, u 3 = ±1, and therefore u\ = u 2 



= 1. 



0. 



u = x 



(l-u 3 )/2 



with a certain 0. In this case, A[i]/y/ci = X (1 ^ us ^ 2 Z(<j)i) 
with a certain angle fa. Therefore, A does not converge to- 
ward £ , since A has eigenoperators with eigenvalues ±1. On 
the other hand, if | r± \ < 1, \\ZUZ'\\ < 1. In this case, 

\\A l \\ < \\UZUZ\\\- 1 /^ 

< \\ZUZ\\ ll/2i 

< ||ZWZ'|| L ' /2J ->0, 

and therefore .4' — > £ in the limit of I — »■ ■. 

This theorem shows that if a resource state is not normal- 
izable even in the limit of infinite size, the MPS matrix is 
always given by A[i]/y/ci = X^- u ^/ 2 Z{<j)i). Apparently, 
such an MPS whose matrix is A[i]/y/ci = X^ 1 "" 3 )/ 2 Z(<j>i) 
cannot allow an arbitrary single-qubit rotation on it. In other 
words, all resource states which are capable of an arbitrary 
single-qubit rotation are at least asymptotically normalizable 
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(i.e. lim A = £). That is, normalizability, which is a suf- 

ficient condition for correct projective measurements in the 
CSs, in the limit of / — > oo is also a necessary condition 
for an arbitrary single-qubit rotation. Accordingly, normal- 
izability (at least asymptotic sense) is a necessary and suffi- 
cient condition for a universal quantum computational wire, 
where one can simulate correct projective measurements and 
arbitrary single qubit rotations. In addition, the theorem com- 
bined with the relationship between the map A and two-point 
correlations (see Appendix B) tells us another interesting fact. 
When lim; A 1 — £ two-point correlations always decay ex- 
ponentially (see Appendix D for details). On the other hand, 
if lim;^ A 1 7^ £, the Kraus operator of the map A is always 
given by A[i]/y/cl = X (1_tl3)/2 Z(0 l ), and hence the two- 
point correlations do not decay. Since any stochastic unitary 
map A can be classified into these two cases due to the theo- 
rem, no MPS of polynomially decaying two-point correlations 
can be described by A of a stochastic unitary map. Thus such 
an MPS cannot be used for MBQC in this framework. 

The above results highlight an interesting relationship be- 
tween computational power and correlation in resource states: 
the many-body states whose two-point correlations do not de- 
cay such as the GHZ and the Ising antiferromagnetic states 
cannot be a universal resource for QCTN. Actually, the MPS 
matrices of the GHZ state are given by A[0] = I/y/2 and 
A[l] = Z/\/2, which apparently does not allow an arbi- 
trary single-qubit rotation in CSs. In addition, a many-body 
state whose two-point correlations exhibit polynomial decay, 
such as the W state, does not allow unitary operations in the 
CS. Accordingly, from the theorem, we can conclude that a 
many-body state whose two-point correlations cannot be up- 
perbounded by an exponentially decaying function, such as 
the GHZ, the Ising antiferromagnetic, and the W states, can- 
not be useful for QCTN, which is the most general framework 
for MBQC to date H33I1 . This also suggests that a ground state 
near the critical point of a quantum phase transition cannot be 
used as a resource state for MBQC, since two-point correla- 
tions of such a state decay as a polynomial function [34]. Evi- 
dence that computational power in MBQC is a robust property 
of a quantum phase has been demonstrated recently on the 
Haldane phase 13511 . Since a two-point correlation is a good 
witness of quantum phase transitions, the present result sug- 
gests that such a robust property should hold in other general 
models of quantum phase transitions. 



V. CONCLUSION 



jective measurements on the finite-size resource states with 
non vanishing two-point correlations by introducing a class 
of MPSs. (iii) We have shown that all resource states which 
are capable of an arbitrary single-qubit rotation are at least 
asymptotically normalizable and exhibit exponentially decay- 
ing two-point correlations. Accordingly, normalizability (at 
least asymptotic sense) is a necessary and sufficient condi- 
tion for a universal quantum computational wire which allows 
simulations of both correct projective measurements and arbi- 
trary single-qubit rotations, (iv) Further, we have seen that a 
many-body state whose two-point correlations cannot be up- 
perbounded by an exponentially decaying function does not 
allow an arbitrary single-qubit rotation, since the matrices of 
such an MPS cannot be described as unitary operators. These 
results (iii) and (iv) have been obtained by showing the theo- 
rem about convergence of a stochastic unitary map. This theo- 
rem itself seems to have an interesting meaning in the context 
of quantum channel. The detailed study is, however, out of 
our scope and an interesting topic for future work. 

Let us finally mention the extension of the above results to 
tensor network states. Since simultaneous single-qubit mea- 
surements are enough to readout results of quantum com- 
putation, we can apply the present argument for tensor net- 
work states by reducing them to MPSs before the readout. 
Thus similar results would be obtained also for tensor network 
states. 

We believe that the present results give us a clue to find 
a novel resource state for QCTN and also help us to build a 
bridge between quantum information science and many-body 
physics, such as the physics of quantum phase transitions. 
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Appendix A: Normalizability and completely depolarizing 
channel 



We have investigated relationships between computational 
power and two-point correlations in the resource states for 
QCTN considering possibilities of projective measurements 
and arbitrary single-qubit rotations, both of which are impor- 
tant properties of a universal resource. Our findings are as 
follows: (i) We have shown that the normalizable resources 
with a finite I exhibit exactly vanishing two-point correlations 
with a finite distance like the cluster state, and always allow 
correct projective measurements in the CS although they are 
of finite size, (ii) We have explored how to perform pro- 



If an MPS is normalizable, 



fn-k(\R), \i/>)) = Tr[|i?)(i?L4°("-%)(V>| 



for all \R) and \tp). This means that 



A ^p=-Tr[p] 
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for all p. Then A°^ n fe ) is a completely depolarizing channel, 
because 

A°^p = ^Tr[p] 



i 

= \Y, u \j)(i\VpV^{j\tf 



Appendix C: Matrix representation of CPTP map 

Consider a unital CPTP map Ap = Y2i A[i]pA '[i], where 
j4t[i]^4[i] = I and A[i]j4J[j] oc /. Each Kraus operator 
can be decomposed into the Pauli matrices as 



A[i] y2 n " rr J 



The CPTP condition reads 



jJ^UaiVpV^iU^. 



J2 a 



Appendix B: Normalizability and two-point correlation 

The expectation value of an observable O k O k+r with re- 
spect to \fy(R, L))i can be calculated as 

(OtO b k+ rh 



E a v a ik e jki + ai a* u + a a a* = 0. 

ijk 

The action of A[i] under its conjugation is given by 



(CI) 

(C2) 
(C3) 



= Tr 



O a k O b k+r mR,L)}1{^(R,L)\ 
1 



/ tt (|J2),|i» 



By using this and condition iC2\ , we obtain Aqo — 
Tlij a ij a ij = 1- Furthermore, by using condition dC3t , we 
(R\ A n -WOl+ r A r - 1 0%A k - 1 \L){L\ \R), obtain A ok = £ 2 a lk a* + a l0 a* k + a l3 a*fyi k = 0. Simi- 

larly we can also obtain A k o = 0. 

Next we define a 3 x 3 matrix A by Aij for i,j ^ 0. A 
quantum state p (i.e. an hermitian operator with trace one) is 
represented in terms of the Pauli basis as p = Z/2+X^' y i cr i/2, 
where Vi is real for all i and its purity < 1. Since the action 
of A[i) map an hermitian operator to another hermitian op- 
erator under its conjugation , AjiV.- L is also real for all i and 
j. It requires that A is a real matrix. Since a CPTP map A 
transforms a state p to another physical state Ap without in- 
creasing its purity, ||.Au||/||f || < 1 for all v, which implies 
that 1 1^| | < 1. 



where O k an d O b k+r are arbitrary operators on the fcth and 
(k + r)th qudits respectively, and 

ofp = J^(i\o a k > b \j)A\j] P Ai[i\. 

ij 

Since n > I, the normalization factor becomes 

U\R),\L)) = \- 

When r > I, 

A"'- 1 = A r ~ l ~ l £ 

for a normalizable resource with a finite I. By using this, we 
obtain 

Ii 



(O a k O b k 



2(R\ 
xTr 



A n-{k+r)Qb 



k+r~, 



\R) 



OiA k - 1 \L)(L\ 



Appendix D: Correlation 

Let us consider the correlation 



On the other hand, the expectation value of 0% and O b k+r can between two Hermitian operators O k on site k and O b k+r on 
be obtained respectively as 



(Olh = Tr 
= Tr 



O a k MR,L))?mR,L)\ 
0%A k ~ x \L)(L\ 



and 



Tr 
2(R\ 



oi +r \m,LMim,L)\ 
i 



site k + r (UOglloo < 1, \\O b k+r \\ x < 1), where 

(0>* = (y(R,L)\0\*(R,L)K. 
Let us define the (super)operator 

d-l d-l 
i=0 j=0 



As a result, 



(OfcOfc+r)* - (Ofc)*(Ofe +r )* = 

for all r > I. 



Then 

(o a k ol 



fn(\R),\L)) 



(R\ 



A^^+^Ol+rA 1 *- 1 



O a k A k - l \L)(L\ \R). 



8 



Let us write 



and 



O a k A k - 1 \L)(L\=a(O a k )- + E(O a k ), 

where a{O k ) is a real number and E{O k ) is a trace-less Her- 
mitian operator. Then, 

A r-l a A k-l | L) (£ | = a(0 a } / + ^^(Qa). 

Note that 

||£(O fe °)lloo < WOtA^imiW^ + ^aiO^ 

_ i 



< d E 



< rf 2 + -| a (0£)|, 
where p is a state and we have used 



(*\O a k \j) 



y^o k {i\o k ){o k \j) 
k 

= d. 

(here, Ofc and \ok) are eigenvalues and eigenvectors of O k ) and 
XI Vci y/cj < d - 



This is because 



? .(-^-[l>(4=) 



> o. 



Thus 



q(Og) 
2/„(| J2>, |i» 
1 

+ M\R),\L)) 
Ci 



li?> 



- 2/„(|i?),|L)) / n (|i2),|L))- 
From this, we also obtain 

«(Og) 



^»-(fc+r) 6 +r / 



2/n(|i2),|i)) 

+ - „„* lrxx (i?|U"- fc £;(og) \r) 



fn(\R),\L)) 

a(O a k ) 

2f n (\R),\L))^ f n (\R),\L)) 



C 3 



{O b k+r )^ 

Vn(\R),\L)) 
1 

+ /n(|i?>,|L)) 

c 4 



" 2/„(|i?),|L» f n (\R),\L)Y 
Note that 



A n -( k+r) & k+r A r - 1 E(I) \R) 



/ n (|i2>,|L» = (fl|(^|L)<L|)|i2> 



where F is a trace-less Hermitian operator. Therefore, 

1 



fn(\R),\L)) 



(R\A n F\R) 



< U n F\\ a 

< \\Ur\ 



where A is the 3x3 matrix defined in the previous section. 
Furthermore, 



\a{Ol)\ = |Tr(0^ fe_1 |L)(i|) 

= \Tr(j2(i\O a k \j)A\j]^[{\) 

< d 2 , 

where £ is a density operator. And 



(R\ 



A n-(k+r) b +r j 



\R) 



Tr(\R){R\A n - {k+r) O b k+r l) 



= \Tr(r,O b k+r I 



< 



Y,mO b k+r \j)\-\Tv[A^] v A[j} 



1,3 



< d 2 , 
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where 77 is a density operator. Therefore, 



where 77 is a state. 



«(0£) 



c 4 



2/ n (|i?),|L)) 2/ n (|i?),|L))2/„(|fl>,|L» 



|Cx! 



< 



Vl{\R)AL)) 
d 4 



2(i-ina)iioo) 2 2 

In the similar way, 



I C 2 



< 2d 2 ||^- 1 ^(02)|| O0 

< ikiic^iu 



IC3 



< 



4«-*J5(Og)|| 
< _|| ( ^)n^|| o 



In conclusion, the two-point correlation is upper-bounded 
by an exponentially decaying function. 
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